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Superconductivity has again become a challenge following the discovery of 
unconventional superconductivity. Resistance-free currents have been observed 
in heavy-fermion materials, organic conductors and copper oxides. The discov- 
ery of superconductivity in a single crystal of U Ge2, ZrZni and URhGe revived 
the interest in the coexistence of superconductivity and ferromagnetism. The 
experiments indicate that: i)The superconductivity is confined to the ferromag- 
netic phase. ii)The ferromagnetic order is stable within the superconducting 
phase (neutron scattering experiments). iii)The specific heat anomaly associ- 
ated with the superconductivity in these materials appears to be absent. The 
specific heat depends on the temperature linearly at low temperature. 

I present a review of the recent experimental results and the basic theoret- 
ical ideas concerning ferromagnetic superconductivity (FM-superconductivity) 
induced by ferromagnetic spin fluctuations. A particular attention is paid to 
the magnon exchange mechanism of FM-superconductivity. 

1 Introduction 

The discovery of unconventional superconductivity caused an explosive growth 
of activities in various fields of condensed-matter research, stimulating not 
only studies of the basic mechanisms leading to this phenomenon, but also 
a widespread search for new technological applications. Resistance-free cur- 
rents have been observed in heavy-fermion materials CeC^^^H], C^-Bei3[2], 
t^-P^ISl, and Ui- x Th x Bei3^\, in organic conductors (SJEHZji copper oxidesdE] 
and layered ruthenate S^-RuC^^O!- The great interest in the unconventional 
superconductivity is in particular due to their rather different normal- and 
superconducting-state properties. The mechanism of superconductivity and 
symmetry of the order parameter are the main puzzling of on-going research. 

The Cooper pairing of conducting electrons is characterized by the gap func- 
tion A CTjCr /(p), which is a 2 x 2 matrix- valued function of the wave vector p. 
There is a symmetry relation which follows from the anti-commutation of spin 
i fermions 

A CT ,^(p) = -A CT> (-p), (1) 
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The inversion symmetry is expressed by 



PA CT , CT ,(p) = A CT>CT ,(-p), 



(2) 



and the time reversal symmetry by 



TA ff , CT (p) = A*_ CT ._ CT (-p), where (a) = (|, j) and (-a) = (j, t) 



(3) 



One can represent the gap matrix in the form 

A CT , ff '(p) = iA(p)(r 2 ) CT)CT ' + id /J (p)(r Al T 2 ) CT)CT /. 



(4) 



where ti,T2 and T3 are the Pauli matrices, A(p) is spin singlet function and 
d/i(p) is spin triplet one. The singlet part of the gap is a symmetric function 
A(— p) = A(p), while the triplet part is an antisymmetric one p) = 



In conventional superconductors, the quasi-particles form Cooper pairs in 
a spin-singlet state A(p) 7^ 0, d fl (p) = 0, which has zero total spin. The 
existence of the gap in the quasi-particle spectrum leads to unusual thermo- 
dynamic properties of the systems: i) The specific heat decreases exponen- 
tially as exp (— A(0)/fcsT) at low temperature, as opposed to the linear tem- 
perature dependence in the Fermi liquid theory [TT], ii) The same anoma- 
lous temperature dependence shows the paramagnetic susceptibility in s-type 
superconductors^ . In the case of the singlet s-pairing, the nuclear spin- 
lattice relaxation rate 1/Ti exhibits a peak just below the superconducting 
transition temperature|13| . Finally, the time-reversal and parity symmetries 
are not broken in conventional superconductors. All these properties are well 
understood on the basis of the Bardeen, Cooper, Schrieffer (BCS) theory of 
superconductivity 

Alternatively, the quasi-particles can form Cooper pairs in a spin-triplet state 
A(p) = 0, d M (p) 7^ 0, analogous to the "p-wave" state of paired neutral fermions 
in superfluid 3 ife|15|. At present, the heavy fermion compound UPt^ and lay- 
ered ruthenate S^RuCD^ are the only known spin-triplet superconductors. The 
most direct way to identify the spin state of Cooper pairs is from measurements 
of their spin susceptibility, which can be determined by the Knight shift|16j. and 
from measurements of nuclear spin-lattice relaxation rate 1/Ti, probed by nu- 
clear magnetic resonance and nuclear quadrupole resonance. No change in spin 
susceptibility was observed on passing through the superconducting transition 
in layered perovskite Sr 2 Ru04^^, and Z7Pt3^|E|- The relaxation rate 1/Ti 
measured in Sr 2 Ru04^\ did not show coherence peak. Muon spin-relaxation 
measurements, for the same materials, reveal the spontaneous appearance of 
an internal magnetic field below the transition temperature. The appearance of 
such a field indicates that superconducting state is characterized by the breaking 
of time-reversal symmetry [2D] - The unconventional nature of the superconduc- 
tivity in U Pt% is confirmed by the observation of three superconducting phases 
denoted as A,B, and C in the if (magnetic field)-T(temperature) phase diagram. 
Phases meet each other at a tetracritical point [21] . 



-<^(p). 
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On theoretical ground, the unconventional superconductivity was described 
as a spin-triplet superconductivity: non-unitary spin-triplet superconductivity 
in U Pt$ [22] , and odd pairing state, which is two-dimensional analogue of the 
Balian-Werthamer state of 3 He, in Sr2RuOi j23- The non-unitary spin- triplet 
is defined by spin-1 complex vectorial function d(p) which satisfy d* (p) x d(p) 7^ 
0. A test for odd-triplet pairing was proposed by T.M.Rice and M.Sigrist 23 . 
They consider a sandwich of thin film of Sr2RuO^ between two singlet supercon- 
ductors with higher transition temperature. Above T sc of Sr2RuOt± this system 
should behave like a standard SNS Joscphson junction where the coupling is due 
to proximity-induced singlet pairs in Sr2RuO^. Below T sc , however, the Joseph- 
son coupling should decrease because as triplet pairing appears in S^RuO^ the 
proximity-induced singlet pairing will be suppressed. The anomalous tempera- 
ture dependence of the Josephson effect would confirm the odd-parity symmetry 
of the order parameter in Sr2RuOi. 

The discovery of superconductivity in a single crystal of UGe2^Q, ZrZw? [2"K] 
and U RhGe^ZQ revived the interest in the coexistence of superconductivity and 
ferromagnetism. The experiments indicate that: i)The superconductivity is 
confined to the ferromagnetic phase. ii)The ferromagnetic order is stable within 
the superconducting phase (neutron scattering experiments). iii)The specific 
heat anomaly associated with the superconductivity in these materials appears 
to be absent. The specific heat depends on the temperature linearly at low 
temperature. 

The interplay of superconductivity and magnetism has a long history. The 
possible coexistence of superconductivity and ferromagnetism was considered as 
a theoretical possibility for weak itinerant ferromagnets many years ago|27j. In 
systems like ErRh^Bi [2E1 121] and HoMorSk^^. it was observed that s-wave 
superconductivity gives way to ferromagnetism at intermediate temperatures. 
In all these compounds, superconductivity and magnetism originate from dif- 
ferent part of the electron system. In contrast, in UGe2, ZrZn2, and URhGe 
apparently the same band electrons are subject to the ferromagnetic and super- 
conducting instability. 

I present a review of the recent experimental results and the basic theoret- 
ical ideas concerning ferromagnetic superconductivity (FM-superconductivity) 
induced by ferromagnetic spin fluctuations. The review should be considered as 
a " progress report" in which I have attempted to focus on some basic aspects of 
this rapidly evolving field. Another goal of the review is to provide the reader 
with a simple overview of the experimental situation in ferromagnetic supercon- 
ductors, summarizing the basic agreements and disagreements between theory 
and experiment. For additional literature the reader should consult other review 
articles |3*T]. 

The review is organized as follows. In the next section, a review of the 
experiments concerning UGe2, ZrZri2 and URhGe compounds is given. The- 
oretical models to describe superconductivity induced by spin fluctuations are 
presented in the third section. A particular attention is paid to the magnon 
exchange mechanism of FM-superconductivity. 
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2 Ferromagnetic superconductivity-experimental 
status 



The coexistence of superconductivity and ferromagnetism has been studied for 
many years, but it has only recently been demonstrated to occur experimentally 24 
1251 126) . The most surprising fact is that the superconductivity occurs only in 
the ferromagnetic phase. The second surprise is that the specific heat anomaly 
associated with the superconductivity in these materials appears to be absent. 

Recent experiments on single crystaljSlj indicate that UGe2 has the base- 
centered orthorhombic crystal structure (Cmmm) with zigzag chains of nearest- 
neighbor uranium ions. The structure is shown in Fig.l taken from Ref . |33| . 



Figure 1: The base-centered orthorhombic Cmmm crystal structure of UGei- 

In the uranium compounds known as " heavy- fermion systems" the 5/ elec- 
trons are highly localized and interact with fcrmions from more dclocalizcd lev- 
els, giving rise to fermion excitations characterized by large effective masses. In 
U Ge2 , however, the 5/ electrons are more itinerant than in many heavy- fermion 
systems. Specific heat measurements show that 7 coefficient 7 = C(T)/T is 
about 10 times smaller than in conventional heavy-fermion C/-compounds(7 « 
35mJ/K 2 ^Q), which suggests that these electrons behave more like the 3d elec- 
trons in the traditional itinerant ferromagnets such as Fe, Co and Ni. UGei 
differs from the 3d metals mainly in having a stronger spin orbit interaction 
that leads to an unusually large magnetocrystalline anisotropy with easy mag- 
netization axis along a (shortest crystallographic axis, Fig.l). 

At ambient pressure UGe-i is an itinerant ferromagnet below the Curie tem- 
perature T c — 52K, with low-temperature ordered moment of ji s = IA^lb/U. 
With increasing pressure the system passes through two successive quantum 
phase transition, from ferromagnetism to FM-superconductivity at P ~ 10 kbar, 
and at higher pressure P c ~ 16 kbar to paramagnetism^! EHJ The resulting 
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pressure p temperature T phase diagram for UGe2 is shown in Fig. 2(a) taken 
from Ref.|m. 
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Figure 2: (a) The p— T phase diagram of U Gei . T c is the Curie temperature T s is 
the superconducting temperature. (b) Pressure dependence of /i at T = 2.3.ftT(full 
circles). The moment obtained by extrapolating the data from above H x to zero 
(squares) . (c) Pressure dependence of the fields H x and H m of metamagnetic 
transitions (at which dM/dH has a local maximum) at T = 2.3K 



At the pressure where the superconducting transition temperature is a max- 
imum T sc = 0.8K, the ferromagnetic state is still stable with T c = 32K, and the 
system undergoes a first order metamagnetic transition between two ferromag- 
netic phases FM2 — * FA/1 with different ordered moments [21] • The pressure 
dependence of the moment /i at 2.3K is shown in Fig. 2(b) taken from Ref. 36 . 
(The symbol M is used for the magnetization in an external field). 

The survival of bulk ferromagnetism below T sc has been confirmed directly 
via elastic neutron scattering measurements . The specific heat coefficient 
7 = C/T increases steeply near 11 kbar and retains a large and nearly constant 
value[37|. 

The resistivity measurements reveal the presence of an additional phase 
line that lies entirely within the ferromagnetic phase. It is suggested by a strong 
anomaly seen in the resistivity. The characteristic temperature of this transition, 
T x (p), decreases with pressure and disappears at a pressure p x close to the 



5 



pressure at which the superconductivity is strongest (Fig2(a)). The additional 
phase transition demonstrates itself and through the change in the T dependence 
of ^(T)!!!!)! E|- in Fig. 3, taken from|3BJ, the temperature dependence of the 
ordered magnetic moment, is shown. A clear change in the T dependence of 
fJ.(T) occurs at T x . 




10 20 30 40 50 60 
T(K) 

Figure 3: Temperature dependence of the ordered ferromagnetic moment. 
Curves correspond from top to bottom to the pressures indicated in the top 
right corner of the frame. 

The field dependence of magnetization at 2.3K for different p is shown in 
Fig. 4, taken from 36 . For pressure p > p x a large increase in the magnetization 
is observed at a field H x . It is defined as the field at which dM/dH has a local 
maximum. It is plotted as a function of pressure in Fig. 2(c). For p > p c the 
magnetization undergoes a second increase at low field H m corresponding to the 
paramagnetic— > ferromagnetic transition. 

It is important to stress that the longitudinal uniform susceptibility, given 
by the slope dM/dH at H = 0, is small at a pressure p x close to the pressure at 
which the superconductivity is strongest, and increases rapidly with increasing 
of p, while the superconducting transition temperature decreases. This experi- 
mental observation suggests that the longitudinal spin fluctuations suppress the 
formation of Cooper pairs. 

The anomaly at T x is quite similar to that observed in a uranium. For a 
uranium, there is direct evidence that the anomalies are due to the formation 
of charge density wave (CDW), resulting from nesting at the Fermi surface. 
Although band structure calculations 33 38 indicate that a spin- majority Fermi 
surface sheet could become nested as a function of the magnetic polarization, 
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Figure 4: The field dependence of the easy-axis magnetization at 2.3K. Curves 
correspond from top to bottom to p=0, 6. 5, 9. 0,11-1, 12. 8, 15. 3, 15. 5, 16. 0,16. 7, 17. 3, 
and 18.2 kbar. The inset shows the ac susceptibility. 

no direct evidence for a charge density or spin density wave has been found in 
neutron diffraction studies |3"9"]. 

The studies of poly-crystalline samples of UGe 2 show that T-P phase di- 
agram is very similar to those of single-crystal specimens of t/C^IUIj. This 
result suggests that high-purity specimens with long mean free paths are not 
necessary, at least in the case of UGe 2 , in order to observe superconductivity. 

It is much less known about the ZrZn 2 and URhGe compounds. The fer- 
romagnets ZrZn 2 and URhGe are superconducting at ambient pressure with 
superconducting critical temperatures T sc — 0.29-fT[2Hj and T sc — 0.25-ff [2"6*| 
respectively. ZrZn 2 is ferromagnetic below the Curie temperature T c = 28. 5K 
with low-temperature ordered moment of /i s = 0.17/is per formula unit, while 
for URhGe T c = 9.5K and fi s = 0.42^. The quadratic low-temperature 
dependence of the squared magnetization is characteristic of simple itinerant 
ferromagnetism 41 . The low Curie temperatures and small ordered moments 
indicate that compounds are close to a ferromagnetic quantum critical point. 

The physical properties of URhGe at zero pressure closely resemble those 
of UGe-i at the high pressures where superconductivity is found. Although 
the space groups and detailed structures of URhGe and UGe 2 are different, 
both structures are orthorhombic and contain zigzag chains of nearest-neighbor 
uranium ions. Neutron scattering experiments reveal that the magnetization is 
almost entirely attributable to uranium 5/ electrons. 
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The observation of a large jump in the specific heat, at the same temper- 
ature as the onset of superconductivity, demonstrates that the transition to 
superconductivity in U RhGe is a bulk phase transition. It also shows that su- 
perconductivity involves the same itinerant electrons that are responsible for 
ferromagnetism. At low temperature the specific heat coefficient 7 is twice as 
smaller as in the ferromagnetic phase. 

The superconductivity in ZrZrii has a number of remarkable features. First, 
it only appears to occur in high-purity single-crystal sample. Second, there is 
no superconducting anomaly in the specific heat. It means that the supercon- 
ducting state is strongly gapless with large portion of the Fermi surface, or even 
all of it, surviving in the superconducting state. Third, in contrast to /7-based 
compounds, the bands at the Fermi energy in ZrZni are predominantly ZrAd, 
and magnetism and superconductivity result from the same Ad electrons. 

The observation of superconductivity in UGe2 , URhGe and ZrZni suggests 
that superconductivity could occur more generally in ferromagnets. The coexis- 
tence of superconductivity and ferromagnetism may thus be more common and 
consequently more important then hitherto realized. 



3 Theory of ferromagnetic superconductivity 

3.1 Lattice theory of strongly correlated systems 

It is generally agreed that the origin of ferromagnetism lies in the Coulomb 
interaction between electrons. The Hamiltonian for electrons with spin a in- 
teracting via Coulomb interaction V(r — r') in the presence of ionic potential 
V ton {r) has the formgU 

H = H + H mt (5) 



where 



2m 



Mr) (6) 



H int - / dV / d\'V{v-v')n a {v)n^(v'). (7) 

aa' 

In equations © and Q ^ CT (r), (r) are electron field operators and h a (r) = 
$+(r)4 r cr (r) is the local density operator. The lattice potential V lon (r) leads to 
a splitting of the parabolic dispersion into bands. The non-interacted problem is 
then characterized by the Bloch wave functions $ Q ,p(r) and the band energies 
£ a(p)i where a is the band's index. One may introduce Wannier functions 
localized at the atomic position by the relation 

= 7^E e ~ lRRs<i> «.p( r )> ( 8 ) 

where iV is the number of lattice sites, and to define the creation and annihilation 
operators ia , c at i a for electrons with spin a in the band a at site Ri by the 
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equations 



c«, = J <Prx* a ,i{r)i> a {r) < — >i(r)=^x«iWw. (9) 

ia 

c + ma = J rf 3 r Xa ,i(r)#+(r) «— > *+(r) = ]Tx^(r)c+ ff . (10) 

After some algebra the Hamiltonian JSJ may be rewritten in lattice representa- 
tion, in terms of creation and annihilation operators 

H = ^""^ t a ijC^ ia .C a j a + — ^""^ Vjimn CgigCpjcr'CSna'Cjmcr, (H) 

qzj'(T af3j5 ijmn aa f 



where the matrix elements are given by 



taij = / d A r X * m {v) 



h 2 

A + V ion (r) 

2m 



Xai(r) (12) 



= J d 3 rj d 3 r' X * al (r) X ^(r')V(r-r') XSn (r') X ^ m (r). (13) 

The Hamiltonian (jllfl is too general to be tractable. One may restrict the 
discussion to a single-band model (a = ,9 = 7 = 5 = 1). The advantage of 
the single-band model is its comparative mathematical simplicity. It is simple 
enough to handle in detail, but yet close enough to physical realty to supply with 
useful information, and the obtained effective model to be of general application. 
It is physically motivated to discuss this model if the Fermi surface lies within a 
single conduction band, and if this band is well separated from the other bands 
and the interaction is not too strong 021 • Accounting for the weak overlap 
between neighboring orbitals in tight-binding description one may restrict the 
sums over all sites of a lattice, to the sum over the nearest neighbors, denoted 
by < i,j >. The remaining one-band, nearest-neighbor Hamiltonian has the 
formO E3 El 

H = -t y (ct a c ja + h.c) + Uy^j^hii + V hihj 

<i,j>a i <ij> 

J ( Si • S, + (14) 



<l,3> 



E 



F {^Mtl^l^ + h - C ) + X E ( £ tr B 3« + h.c) (nj_ CT + Uj-v) . 
<i,j> L a 

where hi = h^ + and Stf = 1/2^2 , • Cicr T o-o-' ^iv '■ The local(Hubbard) 
term describes the Coulomb repulsion (U > 0), the J-term corresponds to the 
direct Heisenberg exchange which is generically ferromagnetic (J > 0) in nature, 
the V-term describes the density-density interaction, the X-term is a density 
dependent hopping and, finally, the F-term describes the hopping of local pairs 
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consisting of an up and down electron. Coulomb repulsion U is the largest 
energy scale in the problem. In a final step one may neglect the nearest-neighbor 
interactions to obtain the simplest model of strongly correlated electrons, the 
Hubbard model 



3.2 Magnon-paramagnon effective theory 

Theories of weak ferromagnetic metals have been developed by several theorists 
ESI EHl S3 ESI EH1 EH EOl EU- The spectrum of the spin excitations has been 
found. It consists of spin fluctuations of paramagnon type and a transverse 
spin-wave branch. The present subsection is devoted to the derivation of an 
effective magnon-paramagnon theory starting from a microscopic single-band 
lattice model of ferromagnetic metals with Hamiltonian 

H = -t ^2 (c+ c ja + h.c.) - J ^ Si ■ Sj + U ^flifhii - fi^nj, (15) 

<i,j>,a <i,j> i i 

where fj, is the chemical potential. 

Given that the Hubbard coupling U is the largest energy scale in the theory 
it is desirable to diagonalize the Hubbard term. We should also define spin-wave 
excitations such that they are the true Goldstone modes. To accomplish both 
of these goals one introduces Schwinger-bosons (0i,<?, <i>\o) and slave- fermions 
(hi, h\, di, d\) represetation for the operators 

c iT = h\ip a + (p] 2 di, Cii = h\ipi2 - <p\ x <k, 

4t^t4iAi = d\dii $i = § fitrTw'tPiv' (16) 

hi = l- hfhi + dfdi, 'flakier + <*-di + h\hi = 1 

The partition function can be written as a path integral over the complex 
functions of the Matsubara time r yio-(r) (^^(t)) and Grassmann functions 
hi(r) (hi(T)) anddi(r) (£(r)). 

Z(l3) = J Dfx(<p,ip,h,h,d,d,)e- S . (17) 
The action is given by the expression 



S= dr 



vPi<r(r)<Pi,r{T) +hi(r)hi(T) +di(r)di(T)^ + H (ip,ip,h,h,d,d) 

i 

(18) 

where /3 is the inverse temperature and the Hamiltonian is obtained from Eas. (|15|) 
and (|16|) replacing the operators with the functions. In terms of the Schwingcr 
bosons and slave- fermions the theory is U(l) gauge invariant, and the measure 
includes S functions that enforce the constraint and the gauge-fixing condition 
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D»(<p,<p,h,h,d,<£) = n 2^ n Dhi ^ Ddi ( r ) 

i,T,a ir 

Y[8 (<Pia(T)<Pia(T) + hi(r)hi(r) + (r)d, (r) - l) JJ 6 (g.f) . (19) 

zr ir 

The ferromagnetic order parameter is a vector field M. The transverse 
spin fluctuations (magnons) are described by (M\ + 1M2) and (M\ — 1M2) 
fields, and the longitudinal fluctuations (paramagnons) by M3— < M3 >. 
Alternatively the vector field can be written as a product of its amplitude 
p = \J M\ + M| + M3 and an unit vector n, M = pn. In the ferromagnetic 
phase one sets M3 =< M3 > +ip and in linear (spin- wave) approximation ob- 
tains p =< M3 > +cp. It is evident now, that the fluctuations of the p field, 
p— < M3 > are exactly the paramagnon excitations in a formalism which keeps 
0(3) symmetry manifest. One can write the effective theory in terms of M- 
vector components, or, equivalently, in terms of p and an unit vector n. I use 
the parametrization in terms of unite vector and spin singlet amplitude because 
the unite vector n describes the true Goldstone modes of the order parameter. 
With that end in view I make a change of variables, introducing new Bose fields 
fia(r) (Mr)) Ref.|52| 

fia(r) = tpia(r) [1 - Ht)K{t) - di(r)di{T)] ~ h , 

jUr) = fe(r)[l-M*W-^(^M] 4 . ( 2 °) 
where the new fields satisfy the constraint 

fia(r)fia(r) = 1. (21) 
In terms of the new fields the spin vector has the form 

Sfr) = lj2fMr^,f^(r) [1 - h z (r)h t (r) - d t {r)d t {r)\ (22) 

aa' 

where nf = fiv T a.a>fi<?' is a unit vector and [l - hi(T)hi(r) - dj(r)di(r)] is 

the spin- vector's amplitude. When the lattice site is empty or doubly occupied 
the spin vector is zero. When the lattice site is occupied by one electron the 
unit vector identifies the local orientation. One can consider the first two 
components nn and m-2 as independent, and then = y/l — — nf 2 . In 
the leading order of the fields, the spin vector has the form 

Sn ~ -na, S i2 ~ -n i2 , - - ~ ~ (hihi + didi) . (23) 

The last equation shows that the longitudinal spin fluctuations are associated 
with the collective fields (hihi + didi). To avoid misunderstandings, it is impor- 
tant to point out that the charge-waves are associated with the collective field 
(didi — hihi) (see the representation of the electron number operator l|16Jl). 
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In terms of the new fields the action has the form 



o 

+ d l {r) 



fia(T)fia{T) + hi{r) ( — - f ia {T)f irT {T) ) hi(T) 



fia(r)fia(r) di(r) 



+ H (/, /, h, h, d,d)\, (24) 



where H (/, /, h, h, d, dA is the Hamiltonian 

To formulate a mean-field theory I drop the terms of order equal or higher 
then six in the Hamiltonian and decouple the four-fermion term, by means of 
the Hubbard-Stratanovich transformation, introducing a real, spin-singlet Si(r) 
field, corresponding to the collective field (Jiihi + didi). Now, the action is 
quadratic with respect to the fermions and one can integrate them out. The 
resulting action depends on the spinons fienfia and the real field Si. It has a 
minimum at the point Si — sq, fi„ = f a , and the stationary condition is 



s = < hihi + didi > 



(25) 



Expanding the effective action around the mean field point one obtains the 
effective model in terms of the spinons / ct (t, r), f a {r, r) and paramagnons <^i(r) 
(2(pi(r) = sq — Si(r)) 53 . The first three terms in the expansion have the form 



S eff - S H + S P 



S 



int> 



(26) 



Sjj is the action of the Heisenberg theory of localized spins. In the continuum 
limit it has the form 



dr / d 6 r 



2Af/ CT (r,r)/ CT (r,r) + ^^^(t.^Mv) 



i/=i 



(27) 

In Ea. l27fl . M = |(1 — So), and sq comes from "tadpoles" diagrams with one h 
or d line, and the renormalized exchange coupling constant J r is calculated in 

m 

Sp is the contribution to the effective action of the paramagnon excitations 
du> d 3 p 



S 



P 



2tt (2tt) 3 



ip(jj, p) r + a h bp <p(-u, -p) 



P 



(28) 



where the constants are obtained from the Lindhard functions for h and d 

fermions in the limit when n and — are small. 

p 

Finally, the spinon-paramagnon interaction has the form 





S mt = M 2 \ JdrJ d 3 rcp(r,r) 



^ d v n(r, r) • 8 v u(t, r) 



(29) 
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where the effective magnon-paramagnon coupling is obtained from triangular 
diagrams with two h and one d lines or with two d and one h lines. 

To analyze the effective model, it is more convenient to rewrite it in terms 
of rescaled spinon fields 



Cx = V2Mf a , Cx = 
The new fields satisfy the constraint 

CrCa = 2M, 

and the action of the effective theory has the form 



(30) 



(31) 



>eff 



dr / d A r 



J, 



C CT (r,r)C CT (r,r) + -f ^ d u M(r, r) • d u M(r, r) 



i/=i 



+ -^(r,r) 



J2d v M(T,T)-d u M(T,T) 



+ S v 



where M is the spin vector 



= M 



(32) 



(33) 



and S p is given by Eq. f2£|l . 

It follows from Eq. (|22|l that the dimensionless magnetization of the system, 
per lattice site is defined by the equation, 



< Sf >-- 



< nf > (1- < /iA + didi >) 



(34) 



At zero temperature < >= 1 and using the Ea. (|25(l one obtains that M 
is zero temperature dimensionless magnetization of the system per lattice site, 
M =< Sf >■ The parameter M depends on the microscopic parameters of 
the theory and characterizes the vacuum. If, in the vacuum state, every lattice 
site is occupied by one electron with spin up, then M = | (s = 0), the 
parameters a and b are equal to zero and r = H- . In this case one can integrate 
over the paramagnons and the resulting theory is the spin ^ Heisenberg theory 
of the localized spins. When, in the vacuum state, some of the sites are doubly 
occupied (< >^ 0) or empty (< hihi >^ 0), then M < \, the relevant 
excitations are the spinon and paramagnon excitations and the effective theory 
is a "spin M" Heisenberg theory coupled to paramagnon fluctuations defined 
by Eas. (|31l32l33|l . The system approaches the quantum critical point when 
M — > (so — * 1), and r(M) approaches zero when M — * 0. Hence, the 
parameter r measures the distance from the quantum critical point. In quantum 
paramagnetic phase (M — 0), the spinon excitations disappear from the spin 
spectrum (see Eqs.lj^ESJ) and one obtains Hertz's effective model^ 4 . (One 
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can add a four-paramagnon term , calculating one- loop diagrams with four h or 
d fcrmion lines, but I have dropped it motivated by the Hertz's result.) 

In thermal paramagnetic phase (above Curie temperature) the spectrum 
consists of spin singlet fluctuations of paramagnon type and spin-i spinon fluc- 
tuations. Well above the critical temperature the spinon has a large gap, and the 
physics of ferromagnetic metals is dominated by the paramagnon fluctuations. 
But just above T c the spinon's gap approaches zero |54| . and the contribution 
of the spin-i fluctuations is essential. 

Below the Curie temperature it is convenient to introduce explicitly the 
magnon excitations a(r, r), a(r, r). To this end, I consider the [7(1) gauge in- 
variant theory 13211 and impose the gauge- fixing condition in the form argC,\ = 0. 
Then the constraint (|31() can be solved by means of the complex field a = (2 
and £1 = \J2M — da. For the components of the spin vector M + — Mi + 
iM2, M~ = Mi — iM2, and M3 one obtains the Holstcin- Primakoff represen- 
tation: 



M + = V2M -da a, M~ = d\j2M -da, 

M 3 = M - da (35) 

The kinetic term in the action and the measure are the same as the kinetic term 
and the measure in the theory of a Bose field. The only difference is that the 
complex fields are subject to the condition 5o< 1. 

In the spin-wave theory one approximates \/2M — da and integrates over 
the whole complex plane. Then, the model is simplified and the effective action 
can be written in terms of magnon a, a and paramagnon (p fields 



'e// 



dui d 3 p r , 2N 

2^ (2tt)3 P ^ ' %UJ + pp > a ^^> 

P) (r + ay + bp 2 ^j <p(-u, -p) (36) 
m\ f -pr dui d 3 pi 

~Y J J-i "2^(27)3 ( - Pl 'PzM^'PiM^'PsM^i -^2,pi -p 2 ) 



where 



p = M J r (37) 



is the spin stiffness constant. 

In the spin-wave approximation the transverse components of the spin fields 
are proportional to the magnon fields 

5+(r,r) = V2Ma(r,r), S"(t, r) = V2Ma(r, r) (38) 

and the field if(r, r) is exactly the paramagnon (longitudinal spin fluctuation) 

5 3 (r,r)- < S 3 <p(t,t). (39) 
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Hence, in Gaussian approximation, spin-spin correlation functions have the form 

D tr {u,p) = -^- 2 , D lon 9 M = 1 (40) 

iuj + pp 2 r + a VA + ^ 

where the longitudinal magnetic susceptibility is 

x = D lon 9(0,0) = - (41) 

3.3 Spin-induced four-fermion interaction 

The next step is to consider model which describes fermions interacting with 
their own collective spin fluctuations. The action of the effective spin-fermion 
theory has the form 



S s _ f = I dr I d 3 r 





4(r,r) f-J-A-. /f )<„(r.r: (42) 



+ ^+(r,r)r: ff ,^(r,r)5"(r,r) 



+ S 



eff> 



where the second term describes the spin-fermion interaction, and S £ jj is the 
action of the effective magnon-paramagnon theory l|3t)|) ■ One may represent the 
spin vector S by means of magnons and paramagnons. The partition function 
can be then written as a path integral over the complex functions of the Mat- 
subara time r a(r, r) , a + (r, r) , tp(r, r) (magnons, paramagnons) and Grassmann 
functions c a (r, r ) , c+ (r, r) 

Z((3) - J D»(a+,a,<p,c+,c <r )e~ 8 *-f. (43) 

In the spin- wave approximation Eas. (|38l39|l the effective action S s _j? is quadratic 
with respect to the spin fluctuations, and one may integrate them out using the 
formula for the Gaussian integral 55 . As a result one obtains an effective Fermi 
theory. It is convenient to write the action as a sum Sj = St_ Q + Stj^ of a 
free part 



Sf- = J dT J4^s [c+(r,p)c CT (r,p)+e CT (p)c+(r,p)c CT (r,p),] , (44) 
o 

where ea(p) are the dispersions of spin-up and spin-down fermions in ferromag- 
netic phase. 

. , p 2 JM . . p 2 JM . . 

e " {p) = ^-»-— ' ^ P) =2^-' l+ — ' (45) 
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and four-fermion interaction resulting from the exchange of spin fluctuations 



J2 

s f-mt = -y 

D lon 9( Xl - x 2 
J 2 



d A xi Id x 2 



c\{x 2 )c^{x 2 ) - c^(x 2 )ci(x 2 ) 



(46) 



d 4 xid 4: x 2 c^(x 1 )c^ (x 1 )D tr (xi — x 2 )c^ (x 2 )ci(x 2 ) 



where x — (r, r), D^ on 9 is the paramagnon propagator, and D^ r is the magnon 
Green function Eq. . 

For the purpose of doing analytical calculations it is convenient to approxi- 
mate the four-fermion interaction with the static one. To this end I replace the 
magnon and paramagnon propagators Eq. 1401 by static potentials 



D tr {u,p) 
D lon 9(u;,p) 



v m (j>) = 2 A 

pp 2 

Vp m (p) = 



(47) 



1 



bp 2 



Let us represent the spin anti-parallel composite held CfCj as a sum of 
symmetric and antisymmetric parts. After some algebra one obtains an ef- 
fective four fermion theory which can be written as a sum of four terms. Three 
of them describe the interaction of the components of spin-1 composite helds 
(TTjTI + IT; ID- The fourth term describes the interaction of the spin singlet 
composite fields fj, — J,|. The Hamiltonians of interactions are 



H 



J 



11 



J 



H U = - 



H 



n 

i 

n 



d 3 Pi 
(2tt) 3 

d 3 p l 
(2tt) 3 



cf(Pi) c 1 h ( p 2) c T(P2 -P3) c l(Pi +Pa) V P m{p 3 ) 



(48) 



C t (Pl) C f ( P 2) C i(P2 - P 3 ) C i(Pl + Pa) V P m{p 3 ) 



(49) 



11 + 11 



[n d3pi 

8 J 11 (2tt) 3 



4 (Pl)c+(P 2 ) + c t (Pl) C f( P 2) 



[ c t(P2 ~ P 3 ) c i(Pi + Pa) + c i(P2 ~ P 3 ) c t(Pi + Pa)] V-(P 3 ) (50) 



where 



J 



n 



d 3 Pl 



4 (Pi)c+(p 2 ) - 4(p x )4(p 2 ) 



16 7 J-- 1 - (2tt) 3 

x [ct(P 2 -P 3 ) c i(Pi +Ps) - c i(P2 -P 3 ) c t(Pi + P 3 )]V+(P 3 ), ( 51 ) 



, . 2M 

v-(p) = —- 

PP 



1 



bp 2 



2M 1 
= — + —7^- (52) 
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The spin singlet fields' interaction Ea. (|51l) is repulsive and does not con- 
tribute to the superconductivity [SB]. The spin parallel fields' interactions Eqs. 
I|48I49|I are mediated by the exchange of longitudinal spin fluctuations and the 
resulting state is spin-parallel Cooper pairing. The interaction of the f J. + J.T 
fields Eq. (|50|) is relevant for magnon- induced superconductivity. It has an at- 
tracting part due to exchange of magnons and a repulsive part due to exchange 
of paramagnons. 



3.4 Paramagnon exchange mechanism of 
FM-superconductivity 

The most popular theory of FM-superconductivity is based on the paramagnon 
exchange mechanism |5 71 155] with Hamiltonians of interaction Eas. l|48|l and (|49|l . 
By means of the Hubbard-Stratanovich transformation one introduces |T an d II 
composite fields, the order parameters, and then the fermions can be integrated 
out. The obtained free energy is a function of the composite fields and the 
integral over the composite fields can be performed approximately by means 
of the steepest descend method. To this end one sets the first derivatives of 
the free energy with respect to composite fields equal to zero, these are the 
gap equations. To obtain the critical temperature T sc one considers the finite 
temperature gap equations linearized at critical temperature [SB - ] 



2Z 2 (0) J (2tt) 3 pm ^ ' e*(k) 

where Z(0) is the mass renormalization constant m*/m — Z(0). For nearly 
ferromagnetic system, it scales with magnetization M like Z(0) ~ ln(l/M)|59|. 
In Eq. (jn3I) , £ct(p) are the spin-tr dispersions l|45|) with m — > m* , and (3 SC is 
the inverse critical temperature. To solve the gap equations, one expands the 
gap functions in spherical harmonics Y\_ m and keeps only the I = 1, m = 
component. Then, the superconducting critical temperature can be obtained 
following McMillan approximation 



1.14w(M)exp[-Z(0)/Af], (54) 



where 



x '-^h k '\ i -iw) v ' m[k) - (55) 

kj are the Fermi wavevectors for the spin-up and spin-down Fermi serfaces, and 
N a (0) is the density of states at the spin-cr Fermi surface. Near the quantum 
phase transition (M — 0), diverges like ln(l/M) and cancels the logarithmic 
singularity of the renormalization parameter Z(0). As a result, the M depen- 
dence of the critical temperature is determined by U)(M). The superconducting 
transition temperatures as a function of the exchange interaction parameter I 
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are depicted in Fig. 5. In the ferromagnetic phase / > 1, and the magnetic tran- 
sitions occurs at / = 1. The parameter r = I — 1 measures the distance from 
the quantum critical point (see Eas. l|4UI41|l 'l and scales with magnetization like 
r ~ M 2 [H1|1>1II3031- 




Figure 5: Paramagnon exchange mechanism of superconductivity. The super- 
conducting transition temperature as a function of the exchange interaction 
parameter I. 

In paramagnetic phase the spectrum consists of spin singlet fluctuations 
of paramagnon type. The exchange of paramagnons leads to a four-fermion 
interaction, which is attractive in triplet channel. The rotational symmetry, 
in continual limit, requires all three components (ft, fj, + J,|, J, J,) of the spin-1 
vector to be nonequal to zero. The superconducting transition temperature as 
a function of the exchange interaction parameter / is depicted in Fig. 5. In the 
paramagnetic phase I < 1, and the parameter r = 1 — I measures the distance 
from the quantum critical point. 

The superconducting critical temperature in Fay and Appel theory increases 
when the magnetization decreases and very close to the quantum critical point 
falls down rapidly Fig. 5. It has recently been the subject of controversial debate. 
It is obtained in 60 , by means of a more complete Eliashberg treatment, that 
the transition temperature is nonzero at the critical point. In |61|. however, the 
authors have shown that the reduction of quasiparticle coherence and life-time 
due to spin fluctuations is the pair-breaking process which leads to a rapid re- 
duction of the superconducting critical temperature near the quantum critical 
point. In order to explain the absence of superconductivity in paramagnetic 
phase of UGe2, URhGe and ZrZrii it was accounted for the magnon param- 
agnon interaction and proved that the critical temperature is much higher in 
the ferromagnetic phase than in the paramagnetic one|fi2|. 

Despite of the efforts, the improved theory of paramagnon induced supercon- 
ductivity can not cover the whole variety of properties of FM superconductivity. 
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3.5 Magnon induced FM superconductivity 

Magnon exchange mechanism of superconductivity has been proposed [53] to 
explain in a natural way the fact that the superconductivity in UGe2, ZrZri2 
and U RhGe is confined to the ferromagnetic phase. The order parameter is a 
spin anti-parallel component ]{ + 11 of a spin-1 triplet (f1\ ti + IT; 44) with 
zero spin projection. The effective Hamiltonian of the system is 

H eff =H + H n+u (56) 

where Hq is the Hamiltonian of the free spin up and spin down fermions with 
dispersions Ea. l|45|) and Hamiltonian of interaction l|5UI) . The transverse spin 
fluctuations are pair forming and the longitudinal ones are pair breaking. 

By means of the Hubbard-Stratanovich transformation one introduces com- 
posite field T4 + 4T an d then the fermions can be integrated out. The integral 
over the composite field can be performed approximately by means of the steep- 
est descend method. To this end one sets the derivative of the free energy with 
respect to composite field equal to zero, this is the gap equation. To ensure that 
the fermions which form Cooper pairs are the same as those responsible for spon- 
taneous magnetization, one has to consider the equation for the magnetization 
as well. 

M = - < c+c T - c+c x > (57) 

The system of equations for the gap and for the magnetization determines the 
phase where the superconductivity and the ferromagnetism coexist. 

The system can be written in terms of Bogoliubov excitations, which have 
the following dispersions relations: 

Ei( P ) = -™ VWW 
E 2 (p)=™-y/e*{p) + \A(p)\* (58) 

2 

where A(p) is the gap, and e(p) = — /i. At zero temperature the equations 
take the form 



M 



\Swf [1 - QirE2m (59) 

A( P ) = ^^nP-k)e ( - £2 (k)) 

8 J (2tt) 3 ^(k) + |A(k)| 2 

The gap is an antisymmetric function A(— p) = — A(p), so that the expan- 
sion in terms of spherical harmonics Yi m (Qp) contains only terms with odd I. 
I assume that the component with / = 1 and m = is nonzero and the other 
ones are zero 

A(p) = A lo (pW-^-cos0. (61) 

V 47T 
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Expanding the potential V-(p — k) in terms of Legendre polynomial Pi one 
obtains that only the component with I = 1 contributes the gap equation. The 
potential Vi(p,k) has the form, 



Vi{p,k) 



3M 

P 
3M 



jr 



Ap 2 k 2 



k- x jp 



p — k 



1 

pk 



ft 



p 2 + k 2 , r' + (p + k f 
■ In 



Ap 2 k 2 r' + (p — k) 2 pk 



(62) 



where 3M/p = 3/p , ft = p/2Mb = p /2b > 1 and r' = r/b « 1. A straight- 
forward analysis shows that for a fixed p , the potential is positive when k runs an 
interval around p (ja — A, p + A) , where A is approximately independent on p. In 
order to allow for an explicit analytic solution, I introduce further simplifying as- 
sumptions by neglecting the dependence of Aio(p) onp (Aio(p) = Aio(p/) = A) 
and setting V\ (pf , k) equal to a constant V\ within interval (p/ — A, pf + A) and 
zero elsewhere. The system of equations I|59I60|) is then reduced to the system 



M 



A 



OO 1 

^ J dkk 2 J dt[l-Q(-E 2 (k,t))} 



J 2 Vt 
32ir 2 



Pf+A 



dkk 2 / dtt 



a e(-E 2 (k,t)) 



Pf—A —1 



e 2 (k) + £-t 2 A 2 



(63) 
(64) 



where t = cos f 



3.5.1 Solution which satisfies aMA < JM 



The equation of magnetization (|63H shows that it is convenient to represent the 
gap in the form A = ^|7c(M) JM, where k(M) < 1. Then the equation 

E 2 (p,t) = 0, (65) 

defines the Fermi surfaces, 



"\J p2 f ± m \j j2M2 ~ -t 2 ^ 2 , Pf = V^lrm (66) 

The domain between the Fermi surfaces contributes to the magnetization M 
in Ea. l|63[l . but it is cut out from the domain of integration in the gap equation 
Ea. l|64|l . When the magnetization increases, the domain of integration in the 
gap equation decreases. Near the quantum critical point the size of the gap is 
small, and hence the linearized gap equation can be considered. Then it is easy 
to obtain the critical value of the magnetization Msc HH • 

When the magnetization approaches zero, the domain between the Fermi 
surfaces decreases. One can approximate the equation for magnetization Ea. (|63|l 
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substituting p^ from Eq. l|66l) in the the difference (Py) 2 — (Pf) 2 an d setting 
p^ = Pf elsewhere. Then, in this approximation, the magnetization is linear in 
A, namely 



A = \J~JkM (67) 
where k runs the interval (0, 1), and satisfies the equation 

/ n 87T 2 , . 

nyl — k -\- arcsin/i = (68) 

rapjj 

The Ea. (|68|l has a solution if mpfJ > 16ir. Substituting M from Ea. (|67|l in 
Ea. l|64|l . one arrives at an equation for the gap. This equation can be solved in 
a standard way and the solution is 



167T ApfK 



3 m 

l 



exp 



3 2Att 



2 v ; PVxmpf 
I(k) = J dtt 2 In (l + Vl - kH 2 



(69) 



Eqs (|67I68I69|I are the solution of the system Eas. l|63l64|l near the quantum 
transition to paramagnetism. 




Figure 6: The zero temperature momentum distribution n, for spin up fermions, 
as a function of q = — and t — cos 9. 

One can write the momentum distribution functions rv(p,t) and n^(p,t) of 
the spin-up and spin-down quasiparticles in terms of the distribution functions 
of the Bogoliubov fermions 

n\p,t) = u 2 (p,t)n 1 (p,t)+v 2 (p,t)n 2 (p,t) (70) 
n l (p,t) = M 2 (p,t)(l-n 1 (p,t))+« 2 (p,t)(l-n 2 (p,i)) 
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where u(p,t) and v(p, t) are the coefficients in the Bogoliubov transformation. 
At zero temperature ni(p,t) = 1, ri2(p,t) = 0(— i?2(p, £)), and the Fermi sur- 
faces Eq. (|65|) manifest themselves both in the spin-up and spin-down momentum 
distribution functions. The functions are depicted in Fig. 6 and Fig. 7. 



1 

n 




Figure 7: The zero temperature momentum distribution n, for spin-down 
fermions, as a function of q = — and t = cos 6*. 

The two Fermi surfaces explain the mechanism of Cooper pairing. In the 
ferromagnetic phase rv and have different (majority and minority) Fermi 
surfaces (see Fig. 8 and Fig. 9, t — graphs). 

n 




Figure 8: The zero temperature momentum distribution n, for spin-up fermions, 
as a function of q = — for t = (the gap is zero) and t = ± (the gap is maximal) . 
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Figure 9: The zero temperature momentum distribution n, for spin-down 
fermions, as a function of q = jj for t = (the gap is zero) and t = ± (the gap 
is maximal). 



The spin-up electrons contribute the majority Fermi surface, and spin-down 
electrons contribute the minority Fermi surface. When the value of the momen- 
tum of the emitted or absorbed magnon lies within interval (p/ — A,pf + A) 
the effective potential between spin-up and spin-down electrons is attracting. 
Hence, if the Fermi momenta pj and lie within interval (pf — A,pf + A) the 
interaction between spin-up electrons, which contribute the majority Fermi sur- 
face, and spin-down electrons, which contribute the minority Fermi surface, is 
attracting. As a result, spin-up electrons from majority Fermi surface transfer 
to the minority Fermi surface and form spin anti-parallel Cooper pairs, while 
spin-down electrons from minority Fermi surface transfer to the majority one 
and form spin anti-parallel Cooper pairs too. As a result, the onset of supercon- 
ductivity is accompanied by the appearance of a second Fermi surface in each 
of the spin-up and spin-down momentum distribution functions (see Fig. 8 and 
Fig. 9, t = 1 graphs). 

The existence of the two Fermi surfaces explains the linear dependence of 
the specific heat at low temperatures: 

C 9ir 2 

T ^ —( N+ W + N-(0)) (71) 

Here 7V ± (0) arc the density of states on the Fermi surfaces. One can rewrite 
the 7 = j: constant in terms of Elliptic Integral of the second kind E(a, x) 

i 1 2s i 1 2s 

(1 + s) 2 E(- arcsin k, -) + (1 - s) 2 E(- arcsin k, -) , 

2i s ~t~ 1 2i s 1 

(72) 

where s = JMm/pl < 1 and k = ^/S/ttA/JM. 



7 = 



mpf 
3k 
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3.5.2 Solution which satisfies \/-A > JM 

In the present sub-chapter one looks for a solution of the system which satisfies 



A > JM 



(73) 



The inequality Eq. (|73l shows that the gap can not be arbitrarily small when the 
magnetization is finite. Hence the system undergoes the quantum phase transi- 
tion from ferromagnetism to FM-superconductivity with a jump. Approaching 
the quantum critical point from the ferromagnetic side, one sets the gap equal to 
zero in the equation for the magnetization (|63|) and considers the gap equation 
(|64() with magnetization as a parameter. It is more convenient to consider the 
free energy as a function of the gap for the different values of the parameter M . 
To this purpose I introduce the dimensionless "gap" x and the parameters s,A 
and g 



3 m 

2' 
*Pf 



Pj 



JM, A 



A 

Pf' 



<J 



J z Vimpf 
8^ 



(74) 



Then the free energy is a function of x and depends on the parameters s, A and 
g. The dimensionless free energy F(x) is depicted in Fig. 10 for A = 0.08, g = 20 
and three values of the parameter s, s = 0.8, s = 0.69 and s cr = 0.595. 



l + A 



F{x) = ^ {T{x) - ^(0)) = x 2 + g J dqq z J dt X 

(s - VO? 2 - 1) 2 + < 2 z 2 ) e( v / (g 2 - 1) 2 + t 2 x 2 - s)- 
(s - V^ 2 - 1) 2 ) e(V(<z 2 -i) 2 - s) 



(75) 



F 




0.5 1 1.5 



Figure 10: The dimensionless free energy F(x) as a function of dimensionless 
gap x. A = 0.08, g = 20, si = 0.8(upper line), S2 = 0.69(middle line) and 
s cr = 0.595(lower line). 
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As the graph shows, for some values of the microscopic parameters A and g, 
and decreasing the parameter s (the magnetization), the system passes trough 
a first order quantum phase transition. The critical values s cr and x cr satisfy 
x cr /s cr = \/3/ft A cr / JM cr > 1 in agreement with Ea. (|73|) . 

Varying the microscopic parameters beyond the critical values, one has to 
solve the system of equations (|63I64H . One represents again the gap in the form 

A = J^k(M)JM (76) 

but now k(M) > 1. Then the equation E2(k,t) = 0, which defines the Fermi 
surface, has no solution if — 1 < t < —1/k(M) and 1/k(M) < t < 1, and has 
two solutions 



pf = ^p 2 f ±m^J 2 M 2 -h 2 A 2 ( 77 ) 

when -1/k(M) < t < l/n(M). 

The solutions (|77|l determine the two pieces of the Fermi surface. They stick 
together at t = ± l/n(M), so that the Fermi surface is simple connected. The 
domain between pieces contributes to the magnetization M in Ea. (|63|l . but it 
is cut out from the domain of integration in the gap equation Ea. H64l) . The 
Fermi surface manifests itself both in the spin-up and spin-down momentum 
distribution functions. The functions are depicted in Fig. 11 and Fig. 12. 




Figure 11: The zero temperature momentum distribution n, for spin up 
fermions, as a function of q = — and t = cos 9. 

When the magnetization approaches zero, one can approximate the equation 
for magnetization Ea. H63|) substituting p^ from Ea. (|77() in the the difference 
(P/ ) 2 — {pj ) 2 and setting p^ = pf elsewhere. Then, in this approximation, the 
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magnetization is linear in A, namely 



A = J-JkM 



(78) 



where k — is the small magnetization limit of k{M). The Ea. (|78|l is a 

solution if mp/J > 16ir (see Eq.$F$). Substituting M from Eq.O in Ea. ljolJ) . 
one arrives at an equation for the gap. This equation can be solved in a standard 
way and the solution is 



16-7T PfA 



exp 



24^ 2 



7T 

4k 3 



(79) 



mp f J 2 V\ 

Eqs i|78!79|l are the solution of the system near the quantum transition to para- 
magnetism. The second derivative of the free energy Ea. (|75|) with respect to 
the gap is positive when mpf,J/16n > (217I-/16) 1 / 3 , hence the state where the 
superconductivity and the ferromagnetism coexist is stable. 




Figure 12: The zero temperature momentum distribution n, for spin-down 
fcrmions, as a function of q — — and t = cos 9. 



The existence of the Fermi surface explains the linear dependence of the 
specific heat at low temperature: 

C 2tt 2 , N 

f = —N(0) (80) 

Here iV(0) is the density of states on the Fermi surface. One can rewrite the 
7 = C/T constant in terms of Elliptic Integral of the second kind E(a, x) 

2s , , 1 „,7T 2s 
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mpf 
3k(M) 



(1 + «)»£?(- 



1 



4 s- 1 



(81) 
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where s < 1. Ea. l81|) shows that for k{M) just above one the specific heat 
constant 7 is smaller in ferromagnetic phase, while for k(M) » 1 it is smaller 
in FM-superconducting phase. The result closely matches the experiments with 
ZrZri2 and URhGe respectively. 

The solutions Eas. (|67l76|l show that magnetization and superconductivity 
disappear simultaneously. It results from the equation of magnetization, which 
in turn is added to ensure that the fermions which form Cooper pairs are the 
same as those responsible for spontaneous magnetization. Hence, the funda- 
mental assumption that superconductivity and ferromagnetism are caused by 
the same electrons leads to the experimentally observable fact that the quantum 
phase transition is a transition to paramagnetic phase without superconductiv- 
ity. 

An important experimental fact is that ZrZn-z and URhGe are supercon- 
ductors at ambient pressure as opposed to the existence of a quantum phase 
transition in U Ge2 ■ To comprehend this difference one considers the potential 
l|t)2|) . The quantum phase transition results from the existence of a momentum 
cutoff A, above which the potential is repulsive. In turn, the cutoff excistence 
follows from the relation [3 = pi 1Mb > 1, which is true when the spin- wave 
approximation expression for the spin stiffness constant p = Mp is used. The 
spin wave approximation correctly describes systems with a large magnetization, 
for example Z7Ge 2 . But in order to study systems with small magnetization, 
one has to account for the magnon-magnon interaction which changes the small 
magnetization asymptotic of p, p — M 1+a po, where a > 0. Then for a small 
M j3 < 1, and the potential is attractive for all momenta. Hence for systems 
which, at ambient pressure, are close to quantum critical point, as ZrZni and 
URhGe, the magnon self-interaction renormalizes the spin fluctuations param- 
eters so that the magnons dominate the pair formation and quantum phase 
transition can not be observed. But if one applies an external magnetic field, 
the magnon opens a gap proportional to the magnetic field. Increasing the 
magnetic field the paramagnon domination leads to first order quantum phase 
transition. 

The proposed model of ferromagnetic superconductivity differs from the 
models discussed in ]§7\ EH] in many aspects. First, the superconductivity is 
due to the exchange of magnons, and the model describes in an unified way 
the superconductivity in UGei, ZrZni and URhGe. Second, the paramagnons 
have pair- breaking effect. So, the understanding the mechanism of paramagnon 
suppression is crucial in the search for the ferromagnetic superconductivity with 
higher critical temperature. For example, one can build such a bilayer compound 
that the spins in the two layers are oriented in two non-collinear directions, and 
the net ferromagnetic moment is nonzero. The paramagnon in this phase is to- 
tally suppressed and the low lying excitations consist of magnons and additional 
spin wave modes with linear dispersion e(p) ~ p|65|. If the new spin-waves are 
pair breaking, their effect is weaker than those of the paramagnons, and hence 
the superconducting critical temperature should be higher. Third, the order 
parameter is a spin antiparallel component of a spin triplet with zero spin pro- 
jection. The existence of two Fermi surfaces in each of the spin-up and spin-down 



27 



momentum distribution functions leads to a linear temperature dependence of 
the specific heat at low temperature. 

The proposed model of magnon-induced superconductivity does not contain 
the relativistic effects, namely spin-orbit coupling which is present in TJGe%. The 
resulting magneto-crystalline anisotropy will modify the spin-wave excitation 
and will add a gap in the magnon spectrum, which changes the potential Eq. 1621) . 

4 Conclusions 

The brief review involves a personal choice of topic and emphasis. Here I have 
concentrated on the coexistence of ferromagnetism and superconductivity. The 
main experimental conclusion is that the same band electrons are subject to the 
ferromagnetic and superconducting instability. The review is devoted to spin 
exchange mechanism of superconductivity. I have not discussed the phonon 
mechanism of triple p-type superconductivity 66, 67 , the symmetry and nodal 
structure of the order parameter |68l IrJ^] and other. 

The physics of ferromagnetic superconductivity involves a subtle interplay 
between magnetism and superconductivity. The concept of magnon induced 
superconductivity is set within the general scheme of itinerant magnetism. It is 
suggested by the fact that superconducting phase lies entirely in the ferromag- 
netic one. The theory explains the ferromagnetic to FM-superconducting quan- 
tum phase transition and the absence of specific heat anomaly. A point which 
requires further theoretical investigation is the presence of an additional phase 
line within ferromagnetic phase. There is no an adequate theoretical explana- 
tion of the resistivity anomaly near the characteristic temperature T x . Further 
theoretical work is needed to understand why T x decreases with pressure and 
disappears at a pressure p x close to the pressure at which the superconductivity 
is strongest. 
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